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Although alterations in gut microbiota composition during acute colitis have been repeatedly 
observed, associated functional changes and the recovery from dysbiosis received little attention. 
In this study, we investigated structure and function of the gut microbiota during acute inflammation 
and recovery in a dextran sodium sulfate (DSS)-colitis mouse model using metatranscriptomics, 
bacterial 16S rRNA gene amplicon sequencing and monitoring of selected host markers. Parallel to 
an increase of host markers of inflammation during acute colitis, we observed relative abundance 
shifts and alterations in phylotype composition of the dominant bacterial orders Clostridiales 
and Bacteroidales, and an increase of the low abundant Enterobacteriales, Deferribacterales, 
Verrucomicrobiales and Erysipelotrichales. During recovery, the microbiota began to resume, but 
did not reach its original composition until the end of the experiment. Microbial gene expression was 
more resilient to disturbance, with pre-perturbation-type transcript profiles appearing quickly after 
acute colitis. The decrease of Clostridiales during inflammation correlated with a reduction of 
transcripts related to butyrate formation, suggesting a disturbance in host-microbe signalling and 
mucosal nutrient provision. The impact of acute inflammation on the Clostridiales was also 
characterized by a significant downregulation of their flagellin-encoding genes. In contrast, the 
abundance of members of the Bacteroidales increased along with an increase in transcripts related 
to mucin degradation. We propose that acute inflammation triggered a selective reaction of the 
immune system against flagella of commensals and temporarily altered murine microbiota 
composition and functions relevant for the host. Despite changes in specific interactions, the 
host-microbiota homeostasis revealed a remarkable ability for recovery. 
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Introduction 

The mammalian intestine is densely populated by 
billions of bacteria belonging to an estimated 
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500-1000 different species predominantly classified 
to the phyla Firmicutes and Bacteroidetes. Under 
conditions of health, host and microbiota live in a 
mutualistic relationship. The host supplies energy 
sources, and the key roles of the microbiota are the 
breakdown of nutrients, the synthesis of hormones 
and vitamins and colonization resistance to patho- 
gens (Nicholson et al, 2012; Stecher et al, 2013). 
Short chain fatty acids (SCFAs such as butyrate, 
propionate and acetate) synthesized during metabo- 
lism of diet- or host-derived carbohydrates are 
actively absorbed and contribute to host energy 
maintenance (Bergman, 1990). Host-microbiota 
interactions are also important at the immunological 
level (Saleh and Trinchieri, 2011). Members of the 
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microbiota or their metabolites, for example, 
Clostridium clusters IV and XIV can induce regulatory 
T-cells, thereby promoting an anti-inflammatory 
immune response (Atarashi et al., 2011). Butyrate, 
which is predominantly produced by Clostridium 
clusters IV and XIV, is a main energy source of 
colonocytes, impacts cell proliferation and differ- 
entiation and lowers the risk of colitis and colorectal 
cancer (Wang et al, 2006; Ploger et al, 2012). 

The human intestine has the highest immune 
activity in the body (Vickery et al., 2011) maintain- 
ing homeostasis with commensals but also respond- 
ing to the invasion of pathogens. In the colon, this 
homeostasis is guaranteed by the combined action of 
a mucin layer, which acts as a first physical line of 
defense, and a tightly regulated network of innate 
and adaptive immunity. Sensing of microbe-asso- 
ciated molecular patterns activates host-signalling 
pathways via pattern recognition receptors, includ- 
ing various Toll-like receptors (TLRs). TLR signals 
activate the synthesis of antimicrobial peptides, 
maintain tight junctions and trigger cytokine production 
in intestinal epithelial cells (Kumar et al., 2009). 

Disturbances of the intestinal homeostasis caused 
by external factors (e.g., dietary shifts and antibio- 
tics treatment) can lead to an imbalance in micro- 
biota composition referred to as dysbiosis (Honda 
and Littman, 2012). Dysbiosis is also associated with 
inflammatory bowel disease (IBD), a chronic inflam- 
mation that is characterized by phases of severe 
(flare-ups) or no symptoms (remission) (Berry and 
Reinisch, 2013). Although the exact etiology of IBD 
remains unclear, several studies suggest that genetic 
predisposition, environmental factors, an overreacting 
immune system and the intestinal microbiota con- 
tribute to the development of chronic inflammation. 

Microbiota alterations frequently observed in 
chronic colitis are an increase of Enterobacteriaceae, 
shifts in abundance of Bacteroidetes and Firmicutes, 
and within-group alterations that can extend to the 
family level or below (Frank et al., 2007; Berry et al., 
2012). Recently, studies applying metagenomics or 
metaproteomics also investigated functional varia- 
tions and concordantly reported alterations in 
metabolic pathways that might affect host- 
microbiota interactions (Erickson et al., 2012; 
Morgan et al., 2012). However, investigations on 
human subjects are mostly limited to the acute 
phase of inflammation. 

With this study we aimed to obtain a more holistic 
picture of the functional links between host and 
microbiota during the time course of colitis devel- 
opment and recovery and to investigate microbial 
functions affected by the inflammation event. We 
used dextran sodium sulfate (DSS) to induce a 
reversible form of colitis (Okayasu et al., 1990). 
Metatranscriptomics was applied to concurrently 
investigate the active microbial community and its 
functional properties, and to link functions to 
bacterial taxa (Berry et al., 2012; Helbling et al., 
2012). To investigate abundance shifts of species- 



level phylotypes, 16S rRNA gene amplicon sequen- 
cing was used. Pro- and anti-inflammatory cytokines 
were monitored as an indicator of the state of the 
immune system. Finally, correlation networks were 
generated applying metatranscriptome, amplicon 
sequence and host inflammation marker data in 
order to identify associations between host and 
microbiota. 

Materials and methods 

Animal experiments 

Wild-type C57BL/6 mice originating from three 
litters and 6-8 weeks of age were provided with 
2% DSS (molecular weight: 36-50 kDa, MP 
Biomedicals, Eschwege, Germany) in autoclaved 
drinking water ad libitum for 5 days, after which 
they received DSS-free drinking water. The DSS 
dose chosen caused a moderate to severe colitis. 
Three mice, which received water throughout the 
time course of the experiment, were used as controls 
for weight determination. Each sampling day, 
representatives from at least two litters were killed. 
The cecum and colon were removed and flushed 
with phosphate-buffered saline to collect luminal 
contents (Berry et al., 2012). Flushed lumen con- 
tents were homogenized, collected by centrifugation 
and either snap-frozen for nucleic acid isolation or 
fixed in 2% paraformaldehyde overnight at 4 °C and 
then stored in 70% ethanol/30% phosphate-buf- 
fered saline at — 20 °C until analysis. The intestine 
was fixed in 2% paraformaldehyde and prepared for 
fluorescence in situ hybridization (FISH), immuno- 
histochemistry and pathology evaluation using 
established methods (Horino et al., 2008). Hematoxylin 
and eosin-stained intestinal colon tissue (photomi- 
crographs taken at 100 x magnification) was scored 
in blind evaluation. For the inflammation score, 
three subscores (inflammation severity, inflamma- 
tion extent and crypt damage) ranging from 0-4 
were given and added together which could give a 
maximum possible score of 12. Animal experiments 
were approved by the institutional ethics committee 
and conducted in accordance with protocols 
approved by the Austrian laws (BMWF-66.006/ 
0002-II/10b/2010). 

DNA and RNA purification and preparation of cDNA 
libraries 

Nucleic acid from lumen contents were extracted 
with a phenol-chloroform bead-beating procedure 
(Griffiths et al., 2000) and separated into DNA and 
RNA using the AllPrep DNA/RNA Mini kit (Qiagen, 
Hilden, Germany), as described before (Berry et al., 
2012). DNA and RNA purified from the same 
extractions were used for 16S rRNA gene amplicon 
sequencing and cDNA generation, respectively. For 
double-stranded cDNA synthesis, total RNA was 
reverse-transcribed using the Superscript Double- 
Stranded cDNA Synthesis Kit (Invitrogen, Lofer, 
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Austria) with modifications as described previously 
(Berry et aL, 2012). Single-stranded cDNA was 
produced from total RNA using Superscript III 
reverse transcriptase according to instructions of 
the supplier (Invitrogen). 



RNA isolation and cDNA synthesis from colon tissue 
For RNA preparation from colon tissue, pieces 
were homogenized in 700 |il RAl buffer of the 
NucleoSpin II RNA isolation kit (Macherey-Nagel, 
Graz, Austria) and processed according to manufac- 
turer's instructions. cDNA was prepared as described 
previously (Stockinger et aL, 2004). 



Metatranscriptome sequencing and data analysis 
Double-stranded cDNA libraries were paired-end 
sequenced using an Illumina HiSeq (Campus 
Science Support Facilities GmbH, Vienna, Austria). 
Read pairs were overlapped using FLASH (Magoc 
and Salzberg, 2011), which yielded between 4.1 and 
15.6 Mio reads of ~170 bp. Metatranscriptomic 
sequencing data were analyzed following an estab- 
lished double RNA analysis pipeline (Urich et aL, 
2008). More than 90% of reads were derived 
from rRNA; the number of putative mRNA reads 
per sample ranged from 75 305 to 1166 930 
(Supplementary Table Si). Community composition 
was determined from 100 000 randomly chosen 
rRNA reads per sample, which were taxonomically 
assigned using CREST (Lanzen et aL, 2012) (bit 
score = 150, top percent = 10, minimal support = 5). 
Putative mRNA tags were compared against the 
NCBI RefSeq database using BlastX and functionally 
and taxonomically classified using MEGAN and the 
SEED functional classification scheme therein (bit 
score = 40, top percent = 10, minimal support = 1) 
(Mitra et aL, 2011). This resulted in 24190 to 
368 950 functionally annotated mRNAs. Metatran- 
scriptomes were generated from three replicates 
collected on days 1, 5 and 14 and from 
four replicates collected on days 8 and 25 
(Supplementary Table S2). PAST (Hammer et aL, 
2001) was used for multivariate analysis of meta- 
transcriptome data. Principal component analysis 
(PCA) was done by eigenvalue decomposition of a 
data correlation matrix. Unpaired Mest and Pearson 
Moment analysis in SigmaPlot 11 (Systat) were 
applied for statistical analysis of variance and 
correlations, respectively. 



16S rRNA gene amplicon pyro sequencing 
Amplicon libraries were obtained using barcoded 
pyrosequencing primers 909F and 1492R, pooled, 
purified and quantified as described before (Berry 
et aL, 2012). Pyrosequencing was performed with 
Titanium reagents on a 454 genome sequencer FLX 
(Roche, Vienna, Austria). All 20 samples collected 
on 5 sampling days were subjected to 16S rRNA 
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amplicon sequencing, but one sample from day 1 
was excluded from analysis due to low sequencing 
depth. Reads were quality-filtered using the ampli- 
con pipeline of the GS Run Processor (Roche) and 
the Pyronoise algorithm in mothur (Schloss et aL, 
2011). Operational taxonomic units were generated 
at 97% sequence identity using UCLUST, as 
described (Berry et aL, 2012). Amplicon sequencing 
libraries, which had a mean number of 17 999 reads 
(range [4415-26 495]) of an average length of 269 nt 
were re-sampled at 3500 reads for a- and (3-diversity 
analyses using QIIME (Caporaso et aL, 2010). 



FISH analysis 

FISH was performed on paraformaldehyde-fixed 
biomass, as previously described (Berry et aL, 
2012). Probes specific for Bacteroides acidifaciens 
(OTU 1812) and Bacteroidales OTU 183 were 
designed based on near full-length 16S rRNA genes 
recovered from screening of clone libraries estab- 
lished with primer pair 8F and 1492R (Lane, 1991). 
Probe specificity was evaluated by checking the 
probes against the SILVA SSU NR108 database 
(Quast et aL, 2013) and stringency of hybridization 
conditions was optimized using either clone-FISH 
(Schramm et aL, 2002) or environmental samples 
containing abundant target populations. For Bacter- 
oides acidifaciens (OTU 1812) probes Bacid87 I 
(S-S-Bac-0087-a-A-24, S'-GCGCCGGTCGCCATCAA 
AAGTTTG-3') and Bacid87 II (S-S-Bac-0087-a-A-22, 
5 / -GCCGGTCGCCATCGGAAGTTTG-3 / ) were used 
with 35% formamide. For Bacteroidales OTU 
183, probe Badal83_437 (S-S-Bada-0437-a-A-22, 
5 / -CGCCCTTTGCTCCCTGACAAAA-3 / ) was used 
with 30% formamide. Probes were deposited at 
probeBase (Loy et aL, 2003). 



Quantification of selected genes and transcripts by qPCR 
Copies of 16S rRNA genes from selected taxa and 
selected microbial transcripts were quantified by 
qPCR using a Mastercycler ep realplex (Eppendorf, 
Vienna, Austria). Reaction mixtures (20|il) contained 
10|il QuantiFast SybrGreen (Qiagen), 1 |il of each of 
the specific primers (final concentration 0.25 jim, 
Supplementary Table S2) and 1 |il of template. 
Running conditions were 95 °C for 5 min, followed 
by 40 cycles of 95 °C for 10 s, annealing for 15 s with 
the temperature specified in Supplementary Table S2, 
and 72 °C for 30 s followed by melting curve analysis. 
Standard curves were generated as described before 
(Berry et aL, 2012). Gene and transcript copy numbers 
were calculated per |ig DNA or RNA, respectively. 
Flagellin expression levels were calculated as the 
ratio of transcript copies and gene copies per |il 
nucleic acid preparation. Unpaired £-test and one-way 
ANOVA combined with an all pairwise multiple 
comparison procedure (Student-Newman-Keuls 
Method) were applied for statistical analysis of 
variance (SigmaPlot 11). 
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For mouse gene expression analysis, qPCRs were 
run with 60 °C annealing temperature on an Eppen- 
dorf realplex cycler in a final volume of 15 |il 
containing 1.5 |il cDNA and 0.3 |iM primer using 
GoTaq qPCR Master Mix (Promega, Mannheim, 
Germany). Primers used are summarized in 
Supplementary Table S3. Significance between 
control (day 1) and test conditions (days 5, 8, 14 
and 25) was calculated using the Relative Expres- 
sion Software Tool (Pfaffl et ah, 2002) software via a 
bootstrapping approach (Pairwise Fixed Reallocation 
Randomisation Test). 2000 random reallocations of 
samples were performed; glyceraldehyde 3-phosphate 
dehydrogenase was employed as a house-keeping 
gene; comparable results were obtained using 
hypoxanthine-guanine phosphoribosyltransferase 
(data not shown). 



Correlation networks 

To construct correlation networks, qPCR data from 
host markers, OTU relative abundances from 16S 
rRNA gene amplicon pyrosequencing and the 
relative abundance of transcripts from SEED cate- 
gories level 1 from metatranscriptome libraries were 
combined. For sequencing-based data, any OTU or 
transcript class that did not reach at least 1% 
relative abundance in any sample was excluded 



from the analysis. Data was log-normalized and the 
Pearson correlation coefficient was generated. Pear- 
son correlation coefficients r>0.6 or r< —0.6 were 
used to construct the correlation networks in R 
software environment, which were then visualized 
using Cytoscape (Shannon et al., 2003). 



Data access 

Amplicon and metatranscriptomic data were 
archived at NCBI Sequence Read Archive under 
SRP008057 and SRX314796, respectively. 

Results 

Host markers of acute colitis 

Development and remission of DSS-induced colitis 
in mice was determined by weight measurements 
and histology scores, and by quantification of 
expression of pro- [IL-lf], IL-6 and IFNy) and anti- 
inflammatory [IL-10) cytokines, the chemokine 
CXCLl and the inducible nitric oxide synthase. 
DSS-treated mice started to lose weight at day 5 and 
reached their lowest weight at days 8 and 9 
(Figure la). At day 14, weight had returned to day 
1 levels, while at day 25 mice had gained 15% 
compared to their starting weight. This weight gain 
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Figure 1 Host parameters during colitis and recovery, (a) Weight curves of untreated (control) and DSS-treated mice. Colitis was 
triggered by the addition of DSS (2%) to drinking water for 5 days (red box); control mice received un-amended drinking water during the 
course of the experiment. Samples were taken before (day 1) and after DSS treatment (day 5), during acute colitis (day 8) and recovery 
(days 14 and 25) (green arrows), (b) Inflammation scores (calculated as the sum of inflammation severity and extend, and crypt depth) 
and area of inflammation (%). (c) Representative photomicrographs (100 x magnification) of hematoxylin and eosin-stained colon 
tissues obtained from each sampling date, (d) Expression ratios of pro- [IL-6, IL1-/3, IFNy and CXCLl) and anti- [IL-10) inflammatory 
cytokines and chemokines and markers of inflammation [iNOS) calculated relative to expression at day 1 and GAPDH. Expression ratios 
were calculated using REST, asterisks (*) indicate significance of expression ratios as determined using Pair Wise Fixed Reallocation 
Randomisation Test (Pfaffl et al, 2002). NS, not significant; ND, not detectable. 
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was comparable to the weight gain of control mice 
not receiving DSS (Figure la). Tissue damage 
peaked at day 8 and was slightly decreased at day 
14; the proportion of affected area declined even 
further to day 25 (Figures lb and c). IL-10 was 
significantly downregulated only at day 5; ILl-fi, 
IL-6 and CXCLl were upregulated (P<0.05) at day 8 
and remained at significantly higher transcriptional 
level at days 14 and 25 (Figure Id). IFNy and 
inducible nitric oxide synthase were significantly 
higher expressed at days 8 and 14 (inducible nitric 
oxide synthase only), but not detectable or not 
significantly different at day 25, indicating resolu- 
tion of acute inflammation (Figure lc). 



Temporal changes in microbiota composition 
16S rRNA transcripts of the metatranscriptomes 
were used to investigate temporal alteration in 
microbial community composition. The bacterial 
phyla Firmicutes and Bacteroidetes contributed the 
majority of 16S rRNA transcripts (together 96-97%); 
with Clostridiales and Bacteroidales being the two 
dominant orders (together 94-95%, Supplementary 
Figure SlA). We applied the 16S rRNA transcripts 
as a measure of relative abundance of bacterial 
groups even though the rRNA content does not 
necessarily reflect cellular abundance and is 
affected in complex ways by the physiological state 
of the cell. PCA of the relative abundance of 
bacterial groups from 16S rRNA transcripts revealed 
that microbiota shifts appeared at the end of DSS 
treatment (day 5) and further alterations occurred at 
the peak of inflammation (day 8) (Figure 2a). At day 
14, the microbial community began to resemble the 
community before inflammation, but did not reach 
its original composition even by day 25. PCA 
furthermore indicated a higher degree of community 
variability between replicates in the inflammation 



and recovery phases (Figure 2a). Shifts from the 
healthy to an inflamed and recovering microbiome 
occurred primarily along the first principal compo- 
nent; the clostridial family Lachnospiraceae was 
associated with the healthy status, while members 
of the Bacteroidetes [Bacteroidaceae, Prevotellaceae, 
Porphyrmonadaceae and Bikenellaceae) were 
indicative of inflammation (Figure 2b). The increase 
of Bacteroidetes 16S rRNA occurred in two stages: 
unclassified Bacteroidales peaked at day 5, while 
members of the Bacteroidaceae, Porphyrmonadaceae 
and Prevotellaceae were most abundant at day 8 
(Supplementary Figure SlC). Some bacterial 
families low in relative 16S rRNA abundance, 
like the Enterobacteriaceae, Verrucomicrobiaceae 
(mainly Akkermansia spp.), Erysipelotrichaceae and 
Deferribacteracae (mainly Mucispirillum spp.), also 
increased during acute inflammation, but 16S rRNA 
reads of these families returned to abundance levels 
of healthy mice at day 14, with the exception of the 
Erysipelotrichaceae (Supplementary Figure SlD), 
which only decreased by day 25. Generally, shifts 
in relative abundance of 16S rRNA of Clostridiales 
and Bacteroidales and the increase of Deferribacte- 
raceae and Enterobacteriaceae 16S rRNA were most 
pronounced in mice that suffered from greatest 
weight loss and the most severe inflammation 
(Supplementary Table S4). During recovery, the 
ratio of Clostridiales to Bacteroidales 16S rRNA 
significantly increased compared to day 1 
(Supplementary Figure SlB). 

In a previous study (Berry et ah, 2012), we showed 
that many microbiota rearrangements during acute 
colitis are hidden when 16S rRNA (gene) sequences 
are grouped on phylum, order or family level. To also 
determine composition shifts of roughly species-level 
phylotypes, which is not possible with rRNA reads 
derived from shotgun metatranscriptomes, 16S rRNA 
gene sequences derived from amplicon pyrosequencing 
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Figure 2 Microbiota dynamics during acute colitis and recovery, (a) Score and (b) loading plot indicating the correlation between 
bacterial communities and the health state of the animals. PCA plots were calculated based on relative abundance of metatranscriptomic 
16S rRNA reads of bacterial families. The score plot shows the relationship of individual samples. The loading plot defines the 
contribution of the original variables (bacterial families) to PCI and PC2 but also indicates the relationship variables to each other, which 
were further determined by Pearson moment correlation analysis, uc, unclassified, 4C0d-2 Cyanobacteria-like lineage. 
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were analysed. Similar to the metatranscriptomic 
data, principal coordinates analysis based on 
weighted and unweighted Unifrac distance indi- 
cated a health-state dependent clustering of bacterial 
communities (Supplementary Figure S2A-D). 
Phylotype richness had already decreased before 
the onset of inflammation (day 5) and only started to 
recover at day 14, whereas community evenness was 
generally not affected (Supplementary Figure S3). 

To investigate compositional dissimilarity within 
bacterial groups across sampling days, Bray-Curtis 
dissimilarity index was calculated based on relative 
abundances or presence and absence of phylotypes 
(Supplementary Figure S4). In general, a shuffling of 
abundances rather than a replacement of commu- 
nity members was the cause for the observed 
community shifts. Shifts were more pronounced 
for members of the Bacteroidetes [Bacteroidales, 
unclassified Bacteroidales and Bacteroidaceae) than 
for the Firmicutes, and appeared already at day 5. 
Bray-Curtis dissimilarity based on relative abun- 
dance of phylotypes assigned to Bacteroidaceae and 
Lachnospiraceae varied to a greater extent at days 
14 and 25 than in healthy mice (day 1). 



Bacterial gene expression during colitis and recovery 
In healthy mice, 69-70% of all putative bacterial 
mRNAs were assigned to Firmicutes (44-47%) and 
Bacteroidetes (23-26%) (53-56% to Clostridiales 
and Bacteroidales). 'Protein Metabolism' and 
'Carbohydrate Metabolism' were the most abundant 
SEED categories; 'Virulence' and 'Motility and 
Chemotaxis' were the third major categories in 
Bacteroidetes and Firmicutes metatranscriptomes, 
respectively (Supplementary Figures S5-7). PCA 
generally clustered mRNA profiles according to the 
day of treatment (Figure 3a). Early changes in gene 
expression were already observed at day 5 but 



became most extensive on day 8 (Figure 3a). mRNA 
reads assigned to the categories 'Motility and 
Chemotaxis' and 'Regulation and Cell Signalling' 
were associated with healthy and recovering 
microbiomes (Figure 3b) and the former category 
correlated positively with Clostridiales and Lach- 
nospiracae (r: 0.516 and 0.528, P<0.05). Associated 
with inflammation were transcripts of 'Virulence', 
'Sulfur Metabolism', 'Stress Response' and 'RNA 
Metabolism' (Figure 3b), all of which positively 
correlated with the Bacteroidales (r: 0.522, r: 0.625, 
r: 0.555 and r: 0.498, respectively, all P<0.05,). 
During acute colitis, transcripts of 'Motility and 
Chemotaxis' were significantly reduced compared 
to day 1, while relative abundances of mRNA 
reads affiliated with 'RNA Metabolism', 'Clustering- 
based subsystems', 'DNA Metabolism', 'Cofactors, 
Vitamins, Prosthetic Groups, Pigments' and 'Fatty 
Acids' were significantly increased (Supplementary 
Figure S5). 



Acute colitis is associated with decreased expression of 
butyrate synthesis genes 

As carbohydrates are major energy sources of the gut 
microbiota, we more closely analysed transcripts 
of the category 'Carbohydrate Metabolism', which 
contributed ~ 20-25% of all transcripts to the 
metatrans crip tome (Supplementary Figure S5). The 
relatively high abundance of transcripts assigned to 
'Fermentation' and its subcategory 'Acetyl-CoA- 
fermentation-to-butyrate' and 'Organic acids' (sub- 
category 'Propionyl-CoA to Succinyl-CoA Module') 
reflected the major fermentative pathways of 
Clostridiales and Bacteroidales, that is, butyrate 
production and succinate-propionate conversion 
(Macy et al, 1978; Pryde et al, 2006), respectively 
(Supplementary Figure S8). Yet, whereas the relative 
abundance of transcripts related to succinate-pro- 
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Figure 3 Shifts in microbiota gene expression during acute colitis and recovery, (a) Score and (b) loading plot indicating the correlation 
between bacterial mRNA expression and the health state of the animals. PCA plots were calculated using relative abundance mRNA 
transcript data of level 1 SEED categories. The score plot shows the clustering of samples obtained on different treatment days. The 
loading plot defines the contribution of the original variables (SEED categories) to the principal components and indicates the 
relationship variables to each other, which was further determined by Pearson moment correlation analysis. CBS, SEED category 
'Clustering-based subsystems'; Cofactors, 'Cofactors, Cofactors, Vitamins, Prosthetic Groups, Pigments'. 
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pionate conversion remained unaffected, those 
associated with butyrate formation were signifi- 
cantly reduced on day 8 (Figure 4, Supplementary 
Figure S8). This decrease occurred within the entire 
metatranscriptome as well as within the metatran- 
scriptome of Firmicutes (Figure 4, Supplementary 
Figure S8). The key enzymes acetyl-CoA 
acetyltransferase, which converts acetyl-CoA to 
acetoacetyl-CoA, 3-hydroxybutyryl-CoA dehydratase 
catalyzing the dehydration of 3-hydroxybutyryl-CoA 
to crotonyl-CoA and butyryl-CoA dehydrogenase 
converting acetacetyl-CoA to butyryl-CoA, were 
significantly reduced on day 8 compared to day 1 
(Supplementary Table S5). Clostridiales possess two 
pathways for the final conversion of butyryl-CoA to 
butyrate (Louis and Flint, 2009). Butyrate kinase and 
phosphate butyryltransf erase, which form butyrate 
from butyryl-CoA via the intermediate butyryl- 
phosphate, were not affected by acute colitis 
(Supplementary Table S5). Transcripts of butyryl- 
CoA:acetate CoA-transferase (ButCoA transferase) 
transfering the CoA of butyryl-CoA to external 
acetate were not detected within the metatranscrip- 
tomes. However, a significant decrease could 
be determined by quantitative PCR (day 1: 
log4.58±0.59, day 8: log3.06 ± 0.35, transcripts 
jig -1 RNA). In confirmation, a decrease of both 
ButCoA transferase transcripts and genes was also 
observed in independent additional experiments 
(Supplementary Figure S9). Additionally, tran- 
scripts associated with the category 'Butanol 
Biosynthesis' were significantly lower on day 8 
(Figure 4). This category was mainly represented by 
pyruvate-formate lyase transcripts of Firmicutes. 
Also significantly reduced was 'Polysaccharide 
Metabolism' due to decreased contribution of 
Bacteroidetes transcripts (Figure 4, Supplementary 
Figure S8). 
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Decrease of clostridial flagellin transcripts and genes 
The SEED category 'Motility and Chemotaxis' con- 
sisted to a large part (80%) of transcripts of clostridial 
'Flagellar Motility'. 'Motility and Chemotaxis' 
dominated the separation of functional clusters 
along the first principal component of the PCA plot 
and was also responsible for the large spread of data 
points during recovery phase (day 25) (Figure 3). 
At days 5 and 8, there were 21 and 49% less 
transcripts assigned to 'Motility and Chemotaxis' 
compared to day 1. (Supplementary Figure S5). 
Likewise, the relative abundance of 'Motility and 
Chemotaxis' was reduced in the metatranscriptomes 
of Firmicutes (Supplementary Figure S7). To deter- 
mine whether this reduction that mirrored the 
decrease of Clostridiales and/or Lachnospiraceae 
16S rRNA was due to species reorganization within 
the phyla Firmicutes or was caused by a down- 
regulation of flagellin gene expression, we additionally 
determined abundances of clostridial flagellin genes 
and transcripts by qPCR (Figure 5). Flagellin genes 
and transcripts were reduced at day 8 compared to 
day 1, and flagellin transcripts were also lower at 
day 25 than at the start of the experiment. Yet, 
whereas at days 1 and 5 the ratio of transcripts to 
genes remained fairly stable, it was significantly 
reduced at days 8 and 25, and there was a strong 
variation in replicates at day 14 (Figure 5) hinting 
that the observed reduction of transcripts during 
acute colitis and at day 25 was not only caused by 
the reduced abundance of Lachnospiraceae but was 
also caused by a downregulation of flagellin 
transcription. 

To further validate the decreased expression of 
clostridial flagellin genes during acute colitis, we 
re-analyzed data from a previous experiment in a 
different genotype (STATl" 7 " mice on C57BL/6 back- 
ground; Berry et ah, 2012), during which we had also 



2.5 




'Carbohydrate Metabolism' 'Fermentation' 



Figure 4 Shifts in 'Carbohydrate metabolism'-related transcripts. Shown are relative abundances of transcripts of the categories 
(a) 'Carbohydrate metabolism' and (b) 'Fermentation' assigned with MEGAN at days 1, 5, 8, 14 and 25 in relation to all mRNA reads. 
Means at days 1 and 8 were compared with unpaired t-test, P<0.05 was considered significantly different (*). 
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observed a decrease of flagellins in metatranscrip- 
tomes. Also here, the ratio of flagellin transcripts to 
genes was reduced ~ fourfold (P<0.05) in mice 
suffering from acute colitis compared to a control 
group (control group [n = 5): 0.09 + 0.18, colitis 
group (i7 = 5): -0.34 + 0.31 [log ratio flagellin 
transcripts/genes]). 



Mucin degradation and utilization during acute colitis 
The SEED category 'Aminosugars', which includes 
(3-hexosaminidases capable of releasing N-acetylgluco- 
samines and N-acetylgalactosamines from O-linked 
mucins, was significantly more abundant in inflamed 




days 

Figure 5 Clostridial flagellin transcript and gene copies during 
colitis and recovery. RNA and DNA were recovered from the same 
nucleic acids preparation, RNA was reverse transcribed into 
single-stranded cDNA. Gene and transcript copies were deter- 
mined with qPCR. The ratio of flagellin to gene copies was 
calculated for individual mice and log-transformed. Means were 
compared by One- Way ANOVA, different letters indicate 
significance (P<0.05). 



mice than in healthy or recovering animals (Figure 4). 
That led us to investigate whether the transcription of 
glycoside hydrolases putatively involved in mucin 
degradation and sulfatases catalyzing the rate-limit- 
ing step of mucin degradation (Wright et al., 2000) 
were affected by inflammation. Beta-mannosidase, 
neuraminidase/sialidase, (3-galactosidase and sulfa- 
tase transcripts were more abundant at day 8 
compared to day 1; oc-mannosidase and (3-hexosami- 
nidase transcripts were significantly enriched 
(Figure 6). Taken together, transcripts of all genes 
presumably involved in mucin degradation were 
significantly higher on day 8 compared to day 1 
(Figure 6). This increase corresponded to a signifi- 
cant increase of transcripts of mucin-degrading 
enzymes within the Bacteroidetes trans criptomes 
during inflammation (day 1: 1.36 + 0.05%, day 5: 
1.39 + 0.26%, day 8: 2.03 + 0.28%, day 14: 
1.31 + 0.44% and day 25: 1.33 + 0.11%). 

Enhanced degradation of mucin could increase the 
availability of mucin components. Indeed, transcripts 
of fucose uptake and utilization were slightly more 
abundant while mannose and N-hexosaminidases 
uptake and utilization remained unaffected 
(Supplementary Table S6). Galactose uptake and 
metabolism was significantly reduced due to the 
decrease of the galactose/methyl galactoside ABC 
transport system MglABC (Supplementary Table S3). 
MglABC is a high-affinity ABC transporter acting as a 
scavenging system in case of low substrate concen- 
tration (Wilson, 1974), thus its downregulation 
suggested increased substrate availability. 



Microbiota-host networks in acute colitis and recovery 
We used correlation networks to identify bacterial 
phylotypes associated with functional transcripts 
and host parameters. IFNy, IL-6, IL-lfi and CXCLl 
expression levels correlated with pathology scores, 
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Figure 6 Mucin degradation-related transcripts during the development and recovery from colitis, (a) Fraction of transcripts related to 
mucin degradation in the metatranscriptomes. (b) Taxonomic assignment of transcripts related to mucin degradation. Means at days 1 
and 8 were compared with unpaired f-test; P<0.05 were considered significantly different (*). uc, unclassified. 



The ISME Journal 



and members of this cluster correlated positively 
with Turicibacteraceae OTU 2651 and negatively 
with Bacteroidales OTU 102 (Figure 7a). Addition- 
ally, inducible nitric oxide synthase expression 
was negatively correlated with Lachnospiraceae 
OTU 2896 and Bacteroidales OTU 1366, and 
IL10 expression was positively correlated with 
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Ruminococcaceae OTU 1695 and Akkermansia 
OTU 711. Examination of the first neighbours of 
phylotypes that correlated with host markers 
revealed extensive intra-phylum discordance in 
phylotype abundances in the Bacteroidetes and 
Firmicutes, for example, Bacteroidales OTU 183 
decreased during inflammation and recovered by 
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Figure 7 Correlation networks of 16S rRNA gene amplicon defined OTUs, functional categories (SEED subsystems level 1) and host 
parameters, (a) Sub-network of inflammation or pathology-associated host markers and their first and second neighbours in microbiota 
data sets, (b) Relative abundance of the inflammation-associated Bacteroides acidifaciens OTU 1812 and the health-associated 
Bacteroidales OTU 183 in amplicon sequencing data and representative fluorescence in situ hybridization images of the two phylotypes. 
(c) Sub-network of first neighbours of TLR5. Black and red edges indicate positive and negative correlations, respectively (r>0.6 or 
r< -0.6). 
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day 25, while Bacteroides acidifaciens OTU 1812 
increased transiently during inflammation 
(Figure 7b), which highlights the inner-group 
dynamics occurring during dysbiosis. As clostridial 
flagellin transcription was heavily impacted during 
inflammation, we also quantified the expression of 
TLR5, the pattern recognition receptor of flagellins 
(data not shown). TLR5 correlated with 'Motility 
and Chemotaxis' and the anti-inflammatory 
cytokine IL-10 (Figure 7c). 

Discussion 

A global perspective on microbial gene expression 
pattern before, during and after acute intestinal 
inflammation has been lacking. We thus performed a 
longitudinal study of the intestinal metatranscrip- 
tome and bacterial community dynamics, and 
selected host parameters in DSS-treated mice to 
address the following questions: How resistant and 
resilient is the gut microbiota to acute intestinal 
inflammation? Which microbial activities are 
affected during acute colitis and recovery? How are 
microbial energy and carbon metabolism impacted 
by the inflammation? What are the impacts of colitis 
on host-microbiota interactions? 

Host and microbiota during progression of DSS- 
triggered acute colitis and recovery 
DSS destroys the mucus biopolymer network result- 
ing in enhanced mucus permeability; however, 



mucus biosynthesis is not affected (Johansson 
et al., 2010). The loss of barrier function allows the 
translocation of microbes, leading to increased 
stimulation of epithelial immune cells and the 
induction of inflammatory responses (Kitajima 
et al., 1999). We started to observe weight loss, 
tissue inflammation and the downregulation of 
IL-10 known to limit the expression of cytokines from 
dendritic cells and macrophages after activation by 
TLRs (Murray, 2005) at day 5 (Figure 8). Microbiota 
alterations occurred concurrently with the onset of 
inflammation. Members of the Bacteroidetes were 
identified as possible indicators of disease onset in 
this study, and were previously shown to be capable 
of inducing colitis in antibiotic-pretreated mice 
(Bloom et al., 2011). However, abundance shifts of 
Firmicutes and Bacteroidetes and intra-phylum 
variations only had a minor impact on microbiota 
gene expression at day 5, highlighting a certain 
degree of functional redundancy of the microbial 
community in the face of slight perturbations. The 
significant upregulation of markers of inflammation, 
extensive weight loss and high inflammation scores 
marked the onset of full-blown colitis (day 8). 

Confirming previous studies, acute colitis was 
characterised by a distinctive dysbiosis (Frank et al., 
2007; Nagalingam et al., 2011; Berry et al., 2012) 
manifested by an increase of low abundant families 
and concurrent intra- and inter-phyla abundance 
shifts of the Bacteroidetes and Firmicutes. 

Acute colitis likewise had major impacts on 
gene expression indicated by temporal abundance 
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Figure 8 Timeline of the interplay of host and microbiota during colitis and recovery. 
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changes of mRNA belonging to broad functional 
categories. With the recovery of the host from acute 
inflammation, severely disturbed microbial traits 
contributing to host-microbiota interactions, such 
as butyrate formation and mucin utilization recov- 
ered quickly to pre-inflammation levels (day 14). 
In contrast, community composition was characterised 
by a memory of disturbance comparable to that 
observed in studies using antibiotics as disturbing 
agents (De La Cochetiere et aL, 2005, Dethlefsen and 
Relman, 2011). Upregulation of pro -inflammatory 
cytokines, inflamed tissue and a compositionally 
altered microbiota persisted up to day 25. In this 
state, the interdependent host-microbiota system 
might remain sensitive towards relapsing flares of 
inflammation. 

Correlation networks were produced in order to 
compare changes in the microbiota with inflamma- 
tion/immune state. This analysis revealed that 
pathology scores and expression of pro-inflamma- 
tory cytokines correlated and also were correlated 
(both positively and negatively) with several bacter- 
ial phylotypes (Figure 7a). With future validation, 
this approach of identifying gut bacterial phylotypes 
that correlate with host health markers may be a 
valuable tool for the development of novel non- 
invasive bacterial biomarkers of intestinal disease 
and may reduce the need for endoscopy. 



Impact of acute colitis on Clostridiales-host 
interactions 

Increasing evidence points to an errant immune 
response directed at commensals during acute 
inflammation (Hand et aL, 2012); one frequently 
identified potential target molecule is the flagellin of 
Clostridiales (Lodes et aL, 2004; Duck et aL, 2007; 
Hand et aL, 2012). A decrease of Clostridiales during 
acute colitis and the concurrent reduction of 
clostridial flagellins as described in this study has 
been repeatedly observed in humans and in animal 
models (Frank et aL, 2007; Erridge et aL, 2010; 
Berry et aL, 2012), but the underlying causes have 
received little attention. Upon disruption of the 
mucus layer and epithelium by DSS treatment, the 
basolateral flagellin receptor TLR5 (Hayashi et aL, 
2001) and intracellular flagellin sensor NLRC4 
(Lightfield et aL, 2008) can be stimulated 
(Figure 8, Rhee et aL, 2005), eventually leading to 
specific antibody production (Vij ay-Kumar et aL, 
2010b). Indeed, we observed a correlation between 
relative abundance of flagellin transcripts and 
expression of TLR5 (Figure 7c). Vijay-Kumar et aL, 
2010a previously showed the regulatory impact of 
TLR5 on the gut microbiota as its inactivation 
resulted in increased bacterial load and changes in 
Firmicutes and Bacteroidetes species composition. 
In colitic mice and human IBD patients, antibodies 
targeting flagellins of commensal strains were 
dominant antigens (Lodes et aL, 2004; Duck et aL, 
2007). This directed immune response against 
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flagellated members of the major intestinal phylum 
Firmicutes might be one reason for the observed 
abundance shift of Clostridiales and Bacteroidetes 
in acute colitis. 

The downregulation of flagellin transcription 
(Figure 5) observed here and in a previous study 
hints that members of Lachnospiraceae have devel- 
oped a strategy to avoid, at least to some extent, 
attacks by the host in a similar way as pathogens that 
tightly regulate flagellin transcription to maintain 
low metabolic costs and to escape recognition by the 
host immune system (Winter et aL, 2010). 
In Bacteroides thetaiotaomicron, the expression of 
a carbohydrate epitope was observed to negatively 
correlate with IgA levels, indicating that commen- 
sals are capable of modifying their surface structure 
in response to host activity (Peterson et aL, 2007). 
Disturbances in flagellin expression could still be 
observed during the recovery phase. This is in line 
with observations that clostridial flagellin-specific 
memory T-cells that were initiated during acute 
colitis caused by T. gondii infection or DSS, 
persisted for more than two months and could 
be rapidly activated in case of another challenge 
(Hand et aL, 2012). 



Impact of acute colitis on Bacteroidetes-host 
interactions 

In healthy mice, the phyla Firmicutes and Bacter- 
oidetes both contributed to transcripts associated 
with enzymes involved in degradation of mucin 
glycoproteins. During inflammation, the Bacteroidetes 
might be the numerically abundant group of 
mucin degraders, as judged from an increase in their 
relative abundance and the expression of the 
respective transcripts. Increased abundance of tran- 
scripts related to mucus degradation was contrasted 
by decreased abundance of 'Polysaccharide Metabo- 
lism', possibly indicating a shift in substrate 
preference within the Bacteroidetes. This shift was 
likely due to shifting intra-phylum composition 
because Bray-Curtis dissimilarity, correlation network 
analysis and FISH revealed that this shift was not a 
common event in all Bacteroidetes phylotypes, but 
rather that a subgroup profited from the inflammation- 
induced changes in the intestinal milieu. One of the 
phylotypes enriched during acute colitis was Bacter- 
oides acidifaciens (OTU 1812, Figure 7b), a species 
originally isolated from mouse cecum (Miyamoto and 
Itoh, 2000). B. acidifaciens was also previously 
identified as an indicator phylotype in acute colits 
(Berry et aL, 2012). B. acidifaciens, together with 
Akkermansia, was recently described as an important 
mucin-degrader in vivo (Berry et aL, 2013). In the 
present study, the increase of B. acidifaciens was 
correlated positively with several mucin glycoprotein- 
degrading enzyme transcripts (oc-l,2-mannosidase, 
r: 0.683, (3-hexosaminidase, r: 0.713, sialidase, 
r: 0.551, P<0.05), suggesting that this strain, as a 
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mucin degrader, might have profited from the changed 
microenvironment during acute colitis. 

In contrast to B. acidifaciens, the metatranscrip- 
tomic analysis did not provide indications for 
Akkermansia as mucus degrader during inflamma- 
tion, although this organism possesses all genomic 
features to degrade complex mucus components 
(van Passel et al., 2011). Possibly, the sequencing 
depth of the metatranscriptome was not sufficient 
for such an identification, since Akkermansia 
abundances were very low in the mice analyzed. 



Temporary alterations of metabolic interactions 
between microbiota and host during acute colitis 
Abundance and intra-taxon shifts of Clostridiales 
and Bacteroidales during the acute phase of inflam- 
mation have far-reaching implications on metabo- 
lite-related host-microbiota interactions. The 
decrease in butyrate-producing Clostridia indicated 
by the decrease of clostridial transcripts of pyruvate- 
formate lyase and key enzymes of butyrate fermen- 
tation pathways, potentially resulted in a reduction 
of butyrate production (Supplementary Figure Si, 
4B, Supplementary Table S4). Concurrently, Bacter- 
oidetes such as the mucolytic B. acidifaciens form 
acetate, succinate and potentially propionate as 
metabolic end products (Macy and Probst, 1979; 
Miyamoto and Itoh, 2000) shifting the balance of 
SCFA towards acetate and propionate. As acetate is 
absorbed and transported to the liver, a shifted 
SCFA balance can lead to a decrease in available 
butyrate as an energy source for colonic epithelial 
cells, but might also affect regulation of colonic cell 
proliferation and differentiation (Wong et al., 2006; 
Ploger et al., 2012). The high abundance of tran- 
scripts of pyruvate-formate lyase, an oxygen-sensi- 
tive enzyme that reversibly converts pyruvate to 
acetyl-CoA and formate in Clostridia (Thauer et al., 
1972) emphasizes its role as a second major 
pyruvate metabolism pathway beside pyruvate: 
ferroxin oxidoreductase (Uyeda and Rabinowitz, 
1971). Pyruvate-formate lyase was also prominent 
in the human intestinal metaproteome (Kolmeder 
et al, 2012). 



Conclusion 

In this study, gut microbiota activity and interac- 
tions with the host during acute inflammation and 
recovery were investigated. With the information 
gained, we proposed a time line of events that 
determined the interplay of host and microbiota 
(Figure 8). Inflammation temporarily had a large 
impact on the microbiota and the host; however, the 
host-microbiota homeostasis revealed a significant 
resilience, particularly evident in the transcriptional 
activity of the gut microbiota. While acute DSS 
colitis is certainly different from IBD, there are some 
similarities such as shift in microbiota composition 



with decrease of flagellated, butyrate producing 
Lachnospiracae and an increase of Enterobacteriacae 
(Frank et al., 2007), as well as changes in mucus 
degrader activity (Png et al., 2010). Therefore, the 
mechanistic outcome of this study might be relevant 
for research on human IBD. 
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